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Abstract 

In the selective withdrawal experiment fluid is withdrawn through a 
tube with its tip suspended a distance S above a two-fluid interface. At 
sufficiently low withdrawal rates, Q, the interface forms a steady state 
hump and only the upper fluid is withdrawn. When Q is increased (or S 
decreased), the interface undergoes a transition so that the lower fluid is 
entrained with the upper one, forming a thin steady-state spout. Near this 
transition the hump curvature becomes very large and displays power-law 
scaling behavior. This scaling allows for steady-state hump profiles at 
different flow rates and tube heights to be scaled onto a single similarity 
proflle. I show that the scaling behavior is independent of the viscosity 
ratio. 



1 Introduction 

A look at Edgerton's photographic sequences of dripping faucets (see, for ex- 



ample, Edgerton et al. (1937)) should instantly convince any skeptic that there 
is something fascinating about a fluid interface changing its topology which 
hints at the richness of the underlying physics. Much attention has been de- 



(see. 


Caflisch & Papanicolaou (1993 


); 


Bertozzi et al. (1994 


); 


Goldstein et al. 


(1993 


); 


Pugh & Shelley (1998)) in the 


same manner as one classifies thermo- 



dynamic transitions. It has been shown that when the topological transition 
involves the formation of a singularity in the fluid flows or interface shapes, as 
is the case with a drop dripping from a faucet, one can use the singularity to 
organize studies of the transit ion systematically (se e, for example, Barenblatt 
(19961) ; ^ister fc Stone (1998| ); [Cohen fc Nagel (200l| )). As the singularity is ap- 



proached, the separation between the length scales characterizing the boundary 
conditions and the length scales characterizing the system grows. Furthermore, 
the fluid systems can enter an asymptotic regime where a subset of the terms in 
the governing equations dominate the flows. The organization of the transition 
studies involves the identification of the different subsets of terms which can bal- 
ance as the singularity is approached. In these asymptotic regimes, the systems 
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typically display scaling of length scales or other physical quantities and often, 
the interfacial profiles near the singularity can all be collapsed onto a universal 
curve. In many systems these singularities manifest themselves in the transition 
dynamics (see , p^jister k, Stone (199^ ); Bensimon et al. (1986); Nagel & Odder 



shede (2000| ); |Zeff et al. (2000| ).) Here, this approach is extended to the study 



of the steady-state interface profiles near the topological transition associated 
with the industrial ly important (see, for example, Muskat (1972); Bear (197^ ; 
Cohen et al. (200l| ); |Ganan-Calvo (199^ )) ])rocess of selective withdrawal. 

In the selective withdrawal experiment a tube is immersed in a filled con- 
tainer so that its tip is suspended a height S above an interface separating two 
immiscible fluids (figure |l]) . When fluid is pumped out through the tube at low 
flow rates, Q, only the upper fluid is withdrawn. The flows deform the interface 
into an axi-symmetric steady-state hump with a stagnation point at the hump 
tip (Fig. 1^). The hump grows in height and curvature as Q increases (or S 
decreases) until the flows undergo a transition where the lower fluid becomes 
entrained in a thin axi-symmetric spout along with the upper fluid. The inter- 
face becomes unbounded in the vertical direction, the stagnation point moves 
from the hump tip into the interior of the lower fluid, and the upper fluid ge- 
ometry becomes toroidal thus changing the topology of the steady state. Once 
the spout has formed, an increase in Q (or decrease in 5*) causes the spout to 
thicken. 

Near the transition, the steady-state mean radius of curvature at the hump 
tip, 1/k, can be orders of magnitude smaller than the length scales characterizing 
the boundary conditions (for example the tube diameter, D). However, for the 
range of parameters explored thus far, even when the system is arbitrarily close 
to the transition, k, while large, remains finite. The manifestation of a true 
singularity would entail the divergence or vanishing of some physical quantities 
or length scales describing the system. In practice, for all physical systems the 
approach to the singularity is cut off at some length scale. For the selective 
withdrawal system where both fluids have a viscosity of about 2 St, Cohen & 



Nagcl (? ) show that the cutoff appears when the mean radius of curvature at 



the hump tip is about 200 /im which is well above the breakdown limit for the 
Navier-Stokes equations. It is still unclear which physical parameters determine 
this relatively large length scale for the hump mean radius of curvature cutoff. 

Nevertheless, fixing S and looking at the steady-state profiles as Q is in- 
creased, Cohen and Nagel observe that, up until the cutoff, both the hump cur- 
vature and height display scaling behavior characteristic of systems approaching 
a singularity. Furthermore, they are able to use the observed scaling relations to 
collapse the hump profiles near the transition onto a universal curve. In doing 
so, they have shown that it is possible to treat this as a "weakly-first-order" tran- 
sition and use the same systematic approach others have successfully applied to 
the study of the drop snap-off problem^. However, understanding how to con- 

^While scaling has been hypothesized by ?^crivos & Lo (1978) for an analogous 3-D system 
it was never shown experimentally. 
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Figure 1: Diagram of the experimental apparatus and photographs of the steady 
state interface at different Q. Fluid is withdrawn from the withdrawal container 
and deposited into a waste container. The upper fluid is then siphoned back 
into the withdrawal container. As shown in the diagram, the parameters for 
this problem include the height of the tube from the interface, S, withdrawal 
rate, Q, surface tension, 7, fluid densities, pu and pi, fluid viscosities, r^u and 
rji, tube diameter, Z), fluid height above the interface, L„, fluid height below 
the interface, Li, container size, and surfactant concentration. Viewed bottom 
to top, the photographs show the evolution of the steady state interface as Q is 
increased. The top photograph of the interface in the spout state is taken from 
Cohen et al. (2001). The middle photograph of the hump at the transition flow 
rate is taken from Cohen & Nagel (?). 
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trol the cutoff curvature is crucial for obtaining a more accurate scaling analysis. 
Moreover, a precise understanding of the transition would allow for control of the 
minimum spout diameter. This control could, in turn be used to advance new 
emerging technologies such as coating micro-particles ([Cohen et al. (2001 )), cre- 
ating mono-dispersed micro-spheres (Ganan-Calvo (1998)), and emulsification 
through tip streaming (Eggers (1997)~ pherwood (1984[ )) which take advantage 
of the selective withdrawal geometry. 

Further insight into this cutoff behavior is gained by comparing with an anal- 
ogous two-dimensional (2- D) problem which roug hly corresponds to replacing 
the tube with a line sink. Jeong fc Moffatt (1992| ) showed that in an idealized 
case where the bottom fluid is inviscid while the top fluid is very viscous, the 2- 
D hump interface forms a two-dimensional cusp singularity beyond sufficiently 
high withdrawal rates. Recently, Eggers (2001 ) showed that when the lower 
fluid has a finite viscosity, the solution to the governing equations changes and 
the system no longer manifests a singularity. Instead, the approach to the sin- 
gularity is cut off and the system undergoes a transition to a different steady 
state. In this new state, a sheet of the lower fluid is entrained along with the 
upper fluid into the line sink. However, the finite lower fluid viscosity prevents 
the hump profiles from scaling onto a similarity solution. In contrast, Cohen| 



& Nagel (?) point out that even though they performed their experiments for 



3-D selective withdrawal with an upper and lower fluid which were comparably 
viscous (about 2 St), they still observe scaling behavior and profile collapse onto 
a similarity solution. Whether or not performing the experiments with a lower 
fluid which is less viscous will allow the system to get closer to the singular- 
ity (as is the case in the 2-D analogue) has so far remained a mystery. In the 
present paper, a systematic analysis of the role which the lower fluid viscosity, 
rji, plays in the selective withdrawal transition shows that neither rji nor the 
viscosity ratio affect the scaling dependence of the hump profiles and the value 
of the cut-off curvature Hu- 



In addition to studying the various scaling dependencies, Cohen fc Nagel 
~(^ ) map out the location of the transition in the S vs. Q parameter space. 
Indeed, this topic is by far the most common feature of the selective withdrawal 
transition discussed in the literature. Yet, even t his seemingly benig n subject 
has its subtleties. For example, the experiments of Blake fc Ivey (1986| ), 



Jirka & 



S. (1979 ) , Harelman et al. (1959| ) , and others used two fluid s which had d ifferent 



densities but were miscible. However, as pointed out by Lister (1989 ), when 
surface tension is absent, there is always some fraction of the lower fluid which 
is extracted so that those experiments were not tracking the actual withdrawal 
transition^. Cohen and Nagel's studies of a pair of immiscible fluids which are 
comparable in their viscosity (about 2 St) are the first to take into account the 
effects of the surface tension in experiments mapping out the transition location 
for low Reynolds number flow. In this paper, these studies are extended to other 
fluid combinations. Once again, I find that neither rji nor the viscosity ratio 

^While these experiments did not track the transition, they may have been sufficient for 
their intended purpose of modeling magma layer mixing during volcanic eruptions. 
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affect the values of S and Q at which the transition occurs. 

The paper is organized as follows: Section || describes the fluids and experi- 
mental apparatus used to make the measurements. This section also describes 
the effects of the container walls and fluid layer thickness on the transition. 
Section ^ addresses the hysteretic nature of the transition. Also addresses is the 
correlation between the amount of hysteresis and the final value for the hump 
curvature, Ku, when the transition takes place at low enough values of S. In sec- 
tion ^ a scaling analysis is performed for a system which is identical to that used 
by Cohen and Nagel with the exception that the viscous lower fluid is replaced 
with one which is two hundred times less viscous. The two sets of results are 
then compared. In section |[ I map out the location of the selective withdrawal 
transition in the S vs. Q parameter space as a function of the upper and lower 
fluid viscosities. Also presented is a comparison of the experimental results with 
the theoretical prediction of Lister (1989 ) for the transition dependence on the 
parameters S and Q. Section |6| presents a discussion of the results. Finally, the 
role of surfactants is addressed in Appendix A. Since it is extremely difficult to 
keep a fluid interface free of surfactants it is important to isolate and control 
their effects. In the current studies the effects on the experimental results are 
found to be rather benign. 



2 Characterization of Fluids and Experimental 
Details 

As shown in figure |^, the parameters important for this experiment are the 
upper and lower fluid viscosities and densities z/;, p„, pi), the interfacial 
tension (7), the oriflce diameter (D), the height of the orifice (S*), the fiow rate 
(Q), the fluid height above the interface (Lu), the fluid layer thickness below the 
interface (ij), the container size, and the surfactant concentration. In order to 
look at the scaling of the steady-state proflles, care must be taken in designing 
an experimental apparatus capable of isolating the proflles near the transition. 
Figure |^ displays a diagram of the apparatus in which the experiments were 
performed. A large tank (30 cm x 30 cm x 30 cm) capable of holding fluid 
layers that were each about 12 cm in height is used as the withdrawal container 
for the fiuids. During the experiments, fiuid is pumped out of the withdrawal 
container and into a waste container. When the system is in the spout state the 
fluids enter the waste container as an emulsion. This emulsion is deposited at the 
bottom of the waste container which is where the droplets comprised of the lower 
fluid remain. The upper fiuid is siphoned back into the withdrawal container 
through large tubes at a rate which matches the withdrawal rate thereby keeping 
constant the upper fiuid layer thickness in the withdrawal container. Note that 
the bottom fiuid layer does not change its thickness when the system is in the 
hump state and decreases its thickness with time when the system is in the 
spout state. However, even for thick spouts (0.1 mm) and for large fiow rates 
(10 ml/sec) the lower fluid layer decreases its thickness at the very slow rate of 
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0.01 mm/min which corresponds to a 0.1% change in the straw height, S. 

The fluids were withdrawn using a B9000 Zenith metering gear pump at- 
tached to a variable speed DC motor. A Dynapar Rotopulser encoder was used 
to read out the withdrawal rate. The pump uses gears to displace fluid from the 
pump intake to the pump outlet. There are small variations in flow rate associ- 
ated with the filling and draining of the gaps between the gear teeth. However, 
at large flow rates, or equivalently at high rotation frequencies, these variations 
damp out and the amplitude of the remaining noise corresponds to a very small 
percentage of the total flow rate. High rotation rates also average out the fluc- 
tuations in the driving motor. Using bigger or smaller pump attachments allows 
for pumping of the fluid at the same flow rate but at a different gear rotation 
rate. This allows for the determination of the effects of noise in the flow rate on 
the transition. A further reduction of the noise in the experiments is achieved 
by siphoning the fluids into the waste container. However, when siphoning, 
the maximum rate of withdrawal (determined by the viscosity, straw diameter, 
transfer tube length, transfer tube diameter, and fluid height difference between 
the withdrawal and waste containers) is substantially smaller. 

The fluids used were heavy mineral oil (HMO), light mineral oil (LMO), sil- 
icone oil (polydimethylsiloxane or PDMS), salt water (Salt H2O), and mixtures 
of glycerin and water (Water/Glycerin). No surface chemistry was observed at 
the two-fluid interfaces even when the liquids remained in contact for periods 
longer than a month. However, a slight change in the transition flow rate at 
fixed S over a period of days indicated that the surfactant concentration at the 
interface was increasing with time. A detailed discussion of the surfactant effects 
appears in Appendix A. The viscosity, v, was measured using calibrated Cannon 
Ubbelohde viscometers immersed in a Cannon constant temperature bath. In 
this manner the viscosity could be determined to within ±5%. Glycerin can be 
diluted with water so that the resultant fluid has 0.01 < f < 9.8 St (Segur & 
Oberstar (1951)). Table |l| lists the values of the viscosities and densities for the 
different fluids. 

The surface tension, 7, of the two-fluid interface was determined using the 
pendant drop m ethod (see, for example, Neumann fc Spelt (1995 ); Hansen & 
Rodsrud (1991 )) which takes advantage of the competition between the surface 
tension and buoyant forces acting on a static drop hanging from a nozzle. The 
buoyant forces distort the drop from a spherical shape. Measuring the distortion 
and density mismatch allows a determination of the surface tension. Implemen- 
tation of this technique on water, toluene, and di-methylformamide showed a 
capability for measuring the surface tensions to within ±10%. Table ^ lists the 
values of 7 for the different fluids. 

The apparatus was illuminated from the rear and a CCD video camera was 
used to image silhouettes of the steady state hump shapes. The images were 
transferred onto a PC where an edge tracing IDL program tracked and recorded 
the points where the derivative of the pixel intensity profile across the hump 
interface was extremized. The profiles were then superimposed onto the original 
images and checked for accuracy. In order to determine the mean curvature at 
the hump tip, k, a Gaussian function was used to fit the tip of the recorded 
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System 2 










System7 


Upper Fluid 


PDMS 


PDMS 


PDMS 


PDMS 


HMO 


HMO 


LMO 


Lower Fluid 


H2O 


Sa.lt H2O 


Wfiter/Glyceriii 


Glycerin 


H2O 


WELter/Glycerin 


H2O 


Su Symbol 


• 


+ 


□ 







T 




7 (dynes/cm) 


43 


40 


29 


23 


35 


31 


34 


P-u (g/ml) 


0.97 


0.97 


0.97 


0.97 


0.87 


0.87 


0.85 


Pi (g/ml) 


1.00 


1.11 


1.24 


1.26 


1.00 


1.24 


1.00 


Ap (g/ml) 


0.03 


0.14 


0.27 


0.29 


0.13 


0.37 


0.15 




10 


10 


10 


10 


1.9 


1.9 


0.57 




0.010 


0.012 


2.3 


9.8 


0.010 


2.3 


0.010 




9.5 XlO~'* 


1.1 XlO~^ 


2.1 XlO~l 


9.3X 10~1 


5.3 X 10~^ 




1.8X 10~^ 




12 ± 3 


27 ± 7 


22 ± 5 


25 ± 5 


31 ± 7 


48 ± 10 


5 ± 1 




0.44 


0.42 


0.41 


0.40 


0.32 


0.30 


0.30 



Table 1: List of the properties for each fluid system studied. Row 1 and 2 
Ust the fluids used. Row 3 hsts the symbols used in plotting the experimentally 
measured Su curves in flgures |, and § Rows 4 through 10 list the values 

of the fluid parameters measured for the different systems. The second to last 
row shows the asymptotic value at which the mean curvature saturates, Husat, 
and the last row lists the power-law exponent a used in fitting the experimentally 
measured Su curves in figures ^ ^ |[ and ^ 

hump profile. The value of k is taken to be the curvature of the fitting function 
at the hump tip. 

Using this apparatus, I tested the effects of the container walls and thickness 
of the upper fiuid layer on the transition flow rate, Qu- Figure |^ shows a plot 
of the transition flow rate Qu vs. the distance Luj from one of the container 
walls while the straw height, S, is held constant. The fluid parameters for this 
particular experiment correspond to those of system 3 in Table |[ When the 
straw is farther than about 2 cm from the container wall, there is no variation 
in the transition flow rate. Figure || shows a plot of the transition flow rate Qu 
as a function of the upper fluid layer thickness L„ for constant S. The fluid 
parameters for this particular experiment correspond to those of system 5 in 
Table |l|. As Lu is increased, Qu first increases but eventually saturates and 
remains constant for Lu greater than about 3cm. In all of the experiments, the 
bottom fluid layer thickness L/ was typically kept at about 12 cm. Measurements 
of the transition flow rates showed no significant variations between systems with 
an Li of 12, 10, and 6 cm. Similar measurements for the Qu dependence on L^,, 
Lu , and Li were conducted for all of the different fluid combinations used in the 
experiments. For the tube diameters {D = 0.16 cm and D = 0.79 cm), tube 
heights (0.07 cm < 5 < 2.0 cm), and flow rates {Q < 20 ml/sec) used in the 
experiments, the container walls were always sufficiently distant and the fluid 
layers were always sufficiently thick so as not to affect the fiows. 



1 



7 



0.20 



0.19 
§ 0.18 



0.17 

3 

0.16 
0.15 



2 4 6 8 10 12 14 
L (cm) 



Figure 2: Plot of the transition flow rate, Qu as a function of the distance 
between the withdrawal tube and the container wall, L^, for system 3 in Table 
|l|. Beyond 2 cm, Qu shows no dependence on Lw Note that the Value of Q„ is 
plotted on a linear scale. 
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Figure 3: Plot of the transition flow rate, Qu as a function of the upper fluid layer 
thickness, L„ for system 5 in Table ^ Beyond 3 cm, Qu shows no dependence 
on Lu. 
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3 Transition Structure and Hysteresis 



In their paper, Cohen fc Nagel ("^ ) map out the location and structure of the 



selective withdrawal transition for a pair of fluids corresponding to system 6 in 
table 11} Here I review these results and note that all of the fluid systems in 
table [^display the same qualitative features. Much of our understanding of the 
selective withdrawal transition can be conveyed by focusing on the effects of 
the parameters S and Q. By flxing Q and tracking the hump proflles the tube 
height below which the interface f orms a spout, Su, can be determined. Figure 
|(taken from [Cohen fc Nagel C^ )) shows that Su oc qo.30±o.05 detailed 
discussion of how this result changes with the various fluid parameters along 
with a comparison to currently available theory is presented in section As 
mentioned in the introduction, for the entire range of parameters explored thus 
far, the evolution of the steady state hump profiles is cut off by the hump to 
spout transition before k diverges. However it is possible that there is more 
than one mechanism that is responsible for the transition cutoff. If this is the 
case it is important to isolate the regimes over which particular mechanisms 
dominate. The cut-off suggests that the transition may be hysteretic. Indeed, 
for transitions occurring at low Q, or cquivalcntly at small Su, the straw height 
at which the spout decays back into a hump is different from Su ■ The difference 
of the two heights or hysteresis is defined as AS*. Figure ^ shows that AS 
decreases exponentially with the ffow rate Q. For this particular system the 
hysteresis can be fit with the function: AS* — 0.04exp^'3/" °^^. For Q > 0.1 
ml/sec, AS was too small to measure. Figure ^ also shows how the hump 
mean radius of curvature, and height, at the transition vary with 

Q. The dramatic decrease in AS* coincides with the onset of a ffat asymptotic 
dependence for at Q > 0.1 ml/sec. In order to quantify this correlation, 

the curvature data is fit with the form = 0.02 + 0.32 exp~'3/° °32 .^^j^j^.]^ 

has the same exponential decay with Q as does the hysteresis. The curvature 
saturation values, Husat, for all of the systems are shown in the second to last 

row in Table |l|. 

A diagram ( iHalc fc Cocak (199l|) ) summarizing the features of the transition 
is shown in figure |5|. The vertical axis labeled q relates how much of the lower 
fluid is withdrawn as a function of the total fluid withdrawal rate Q and straw 
height S. The states accessible to the system are represented by a sheet which 
is embedded in this three dimensional parameter space. The part of the sheet 
which corresponds to the system in the hump state is located in the q = 
plane. An experiment where the flow rate is held fixed and the straw height 
is slowly reduced would correspond to a path which traverses the sheet in a 
direction which is parallel to the S axis. The fold in the sheet represents the 
hysteresis in the transition. As the straw height is reduced, the system runs out 
of hump states and makes a discontinuous jump at S = Su (depicted by the 
dashed arrow pointing upward) to a part of the sheet which has a finite value 
of q. If S is reduced further, the sheet slopes towards higher values of q and 
more of the lower fluid is withdrawn. If on the other hand, S is increased, the 
system remains on the upper branch of the sheet until the spout states run out 
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Q (ml/sec) 



Figure 4: Plots of the transition tube height, Su, hysteresis, AS*, transition 
hump height, hu, and transition radius of curvature 1/ku, as a function of the 
flow rate Q. Su oc q0-3±0 05 (solid line). The AS" data is fit with an exponential 
decay (dashed) of the form: AS" = 0.04exp~'5/o o32 although the range is not 
sufficient to exclude a power-law decay. The mean radius of curvature, 
is fit (dash dot) with the form: 1/k u = 0.02 + 0.32 exp-<3/o n32^ rpj^jg gg^^g 
reproduced from Cohen & Nagel (?). 
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at which point the system makes another discontinuous jump (depicted by the 
downward dashed arrow) to a part of the sheet where q — 0. The separation 
between the two arrows along the S axis corresponds to AS*. As the diagram 
indicates, at large values of Q and S, the amplitude of the fold or equivalently 
AS decreases dramatically. 




Figure 5: Bifurcation diagram summarizing the features of the selective with- 
drawal transition. The vertical axis labeled q depicts the amount of the lower 
fluid withdrawn as a function of the withdrawal rate Q and straw height S. The 
folded two dimensional sheet represents the states which are accessible to the 
system. The fold represents the amount of hysteresis in the transition which 
becomes too small to measure at large enough Q and S. 

It is unclear whether at high enough Q or S the system undergoes a bifur- 
cation where the hysteresis vanishes, or if there always remains some small fold 
in the sheet which is too small to be detected by the experiments. It seems 
unlikely that the transition could have a cut-off (and remain discontinuous) and 
yet have no hysteresis. Such a situation would entail the disappearance of the 
hump solutions at the exact location in the parameter space where the spout 
solutions emerge. It is likely that numerical investigations which are able to 
detect both stable and unstable solutions will be needed in order to answer this 
question definitively. The exponential decay in both the AS and 1/k„ curves 
suggests that even if there does exists a small amount of hysteresis at large S 
or Q, there is a separate and localized mechanism that generates the transition 
hysteresis at low values of S* or Q. 

When transitions occur at low S* or Q, the ratio of the withdrawal tube di- 
ameter, D, over the straw height, 5, is of order one. As a result, the circulation 
rolls set up by withdrawal press against the two-fluid interface and the stream- 
lines near the hump tip become distorted. By changing the flow geometry, the 
straw diameter sets a length scale for the hump mean radius of curvature at the 
transition, 1/ku- Figure ^ shows the effect of a factor of five increase in D on the 
Su and 1/ku curves. For the range of flow rates explored, AS is always about 



11 



an order of magnitude smaller than Su- Therefore, it is not surprising that the 
two Su curves in figure ^ show no significant dependence on D. On the other 
hand, the 1/ku data show that the onset of the flat asymptotic dependence for 
1/ku occurs at higher Q for larger straw diameters. This evidence indicates 
that the hysteresis observed in the experiments results from the finite width of 
the withdrawal tube. At large values of Su these finite size effects vanish. We 
therefore restrict our scaling analysis to this regime. 




10"^ 10"^ 10° 10^ 

Q (ml/sec) 

Figure 6: Dependence of the Su and 1/k vs. Q curves on the withdrawal 
tube diameter, D. The closed symbols reproduce the Su and k„ curves in figure 
^ where D = 0.16 cm. The open symbols depict the values of Su and Ku as 
a function of Q when D = 0.79 cm. The solid line is a fit to the D = 0.16 
data which shows that Su oc Q"-3±o.o5^ rj.^^ 

mean radius of curvatures, 1/ku 
are fit (dash) with the forms: 1/ku = 0.02 + 0.32 exp-'3/0 032 ^j^^ ^ 
0.02 + 2.3exp-<3/o-22 ^j. ^-j^g ^ ^ o.l6 and D = 0.79 cm data sets respectively. 



4 Scaling Analysis And Similarity Solutions 



Although the cut-off in the transition exists even for high values of Q, Coherj 



fe Nagcl (?) showed that the separation in length scales between the hump 



radius of curvature and the boundary conditions is sufficient for the system to 
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display scaling and a similarity solution near the transition. As mentioned in 
the introduction, in the 2-D analogue to the selective withdrawal problem the 
viscosity of the lower fluid, i/;, greatly affected the observed scaling dependencies. 
In order to determine the effects of vi on the transition structure in the present 
3-D problem, a comparison of the scaling behavior for two systems with different 
values of vi is presented below. 

Figure |^ reproduces the Cohen and Nagel scaling curves for the hump cur- 
vature and height for system 6 while figure js] shows the results of an identical 
treatment performed on data taken for system 5 (where vi is two hundred times 
smaller). Figures [T^ and^, show a plot of the mean curvature at the hump tip, 
K (where < k < as a function of Q for six representative data sets each 
of which corresponds to measurements taken at a different fixed straw height S. 
All of the curves display a steep rise in the curvature with increasing Q. The 
insets to figures ^ and ^ show that, as the flow rate approaches Qc (used as a 
fitting parameter), the steep rise in the curvatures takes the form of power-law 
divergences with exponents of about -0.84 and -0.80 respectively. For each fluid 
system, all of the data sets (15 sets for large vi system and 20 for low vi system) 
have the same power-law exponent for the divergence. However, the power-law 
prefactors, Ck(S'), change with S. These prefactors are scaled out in the insets 
for clarity. 

Figures l^b and^D, show a plot of the hump height, h^ax (where < hmax ^ 
hu), as a function of Q for six data sets each of which corresponds to measure- 
ments taken at a different fixed straw height S. The insets to figures |^ and ||b 
show that as the flow rate approaches Qc (obtained from insets to figures ^ and 
||a), the hump heights approache the height he (used as a fitting parameter) as 
power-laws with exponents of about 0.73 and 0.70 respectively. While the curves 
for each fluid system show that hmax approaches h^ with the same power-law 
exponent, the power-law prefactors Ch{S) change with S. These prefactors are 
scaled out in the insets for clarity. 

We combine the two scaling dependencies in figures 0c and ^ and plot 
ha-h^an: ^ functiou of the normalized curvature -. The normalization con- 
stant, n, scales out the power-law prefactors for curves corresponding to dif- 
ferent straw heights. The two figures show that ^'"^""^ scales as (^)^ where 
(} = -0.85 ± 0.09 for the high vi system and (3 = -0^86 ± 0.10 for the low vi 
system. These results indicate that for this range of straw heights, both fluid 
systems show no change in the power-law exponent —fi even though the height 
he and the power-law prefactor n change with S. 

Note that n{S) = Ch{S)[ck{SY]. The constants n{S) decreases roughly as 
exp~^-^'^ for both systems. While both systems have the same form for the de- 
cay of the constants n(S'), in the lower vi system, the values of n{S) are larger 
by about ten percent. Also, as indicated by table |l| , the asymptotic value of 
ku at the high Q regime is different by about a factor of 1.5. A detailed anal- 
ysis of table |l] described in Section 6 suggests that these very slight differences 
are attributable the different values Ap in the two systems. This comparison 
demonstrates that the lower fluid viscosity vi does not affect the scaling relations 
near the selective withdrawal transition. 
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Figure 7: Scaling for the hump mean curvature, k, and height hmax- Figure 
^ plots K, vs. Q for six tube heights. Each curve displays diverging behavior 
with increasing Q. For each tube height, a critical flow rate, Qc, is chosen 
as a fitting parameter so that as Q approaches Qc, the curvatures increase as 
power laws with an exponent of about -0.84 (|^a inset). The prefactors to the 
curvature power laws, Ck(S'), are scaled out in the inset. Figure ^ plots hmax 
vs. Q. For each S, a critical hump height, is chosen as a fitting parameter. 
As the flow rate approaches Qc (obtained from |^a inset), the hump heights 
approach the critical heights as power-laws with an exponent about 0.73. The 
prefactors, Ch{S), to the power laws in the inset are scaled out. Figure |^ 
plots {he — hmax) / hmax VS. k/h for the entire data set corresponding to fifteen 
different straw heights. The prefactor n roughly decreases as exp~^-^"^. The line 
corresponds to a power-law with an exponent of -0.85. This figure is reproduced 



from Cohen & Nagel (?) 
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Figure 8: Scaling for the hump mean curvature, k, and height hmax- Figure^ 
plots the K vs. Q for six tube heights. Each curve displays diverging behavior 
with increasing Q. For each tube height, a critical flow rate, Qc, is chosen 
as a fitting parameter so that as Q approaches Qc, the curvatures increase as 
power laws with an exponent of about -0.80 (||a inset). The prefactors to the 
curvature power laws, c^iS), are scaled out in the inset. Figure ||b plots hmax 
vs. Q for six tube heights. For each S, a critical hump height, he, is chosen as a 
fitting parameter. As the fiow rate approaches Qc (obtained from||a inset), the 
hump heights approach the critical heights as power-laws with an exponent of 
about 0.70. The prefactors, Ch{S), to the power laws in the inset are scaled out. 
Figure^ plots {he — hjnax)/hmax vs. k/ti for the entire data set corresponding 
to twenty different straw heights. The prefactor n roughly decreases as exp"^'^'^. 
The line corresponds to a power-law with an exponent of -0.86. 
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The scaling analysis shows that the hump profiles behave as though they 
are approaching a singular solution where, at the flow rate Qc the hump height 
would be equal to he and the mean curvature, k, would diverge. The transition 
cuts off the evolution of the hump states preventing the system from getting 
arbitrarily close to the singularity and limiting the precision with which the 
power-law exponents can be determined. Qc changes with S indicating that the 
system can approach a continuous line of singularities which closely tracks the 
hump to spout transition curve in the S vs. Q parameter space. 

Cohen & Nagel (7) showed that the scaling relations for the hump height 
and curvature can be used to collapse the hump profiles near the transition onto 
a universal curve . The quantities ^ and ^''=~^'"°°' track how quickly the radial 
and axial length scales decrease as the system approaches the singularity. In 
accordance with the definitions of Cohen and Nagcl, the scaled variables are 
defined as: 

, , hr — h{r) TK , , 

^ h - h ""^^ R^—, (1) 

Here h[r) is the hump profile and the value of he is taken from the scaling 
relations. This transformation shifts and scales the profiles so that the singu- 
larity is located at the origin and the maximum hump heights are located at 
H = 1 and R = 0. Figure ^ reproduces the similarity analysis of Cohen and 
Nagel (system 6) for the hump profiles near the transition while figure |l^ shows 
the results of an identical treatment performed on data taken for the small i^i 
system (system 5). 

Figures ^ and |l^ show a series of eight scaled profiles for the S = 0.830 cm 
data set in the large i^i system and the S = 0.667 cm data set in the small i^i 
system. The bottom insets in the figures show an overlay of the different hump 
profiles which are scaled in the main figures. Both figures display excellent 
collapse for the hump profiles. The solid lines are power-laws which fit the scaled 
profiles in the region beyond the parabolic hump tips. The scaling relation in 
figures and || can be used to predict the exponent, x, in the power-law fits. 
Inserting the observed scaling dependence ^""^"""^ = {^)^ into the observed 
form H{R) = which fits the profile, the following relation is obtained: 

hmax ^ ^ 

Since, for a given r, the functions ^'^ '^^^"^ and have constant values, x must 
equal — /3. More intuitively, as Q is increased, the parabolic tip regions decrease 
their radial length scale and are simultaneously pulled towards the singularity 
in the axial direction leaving behind power-law profiles with exponents that 
reflect the rate at which these length scales are decreasing. The fits to the 
profiles indicate that x — 0.72 ± 0.08 and x — 0.72 ± 0.10 for the high and 
low vi systems respectively. Both of these exponents are within error (although 
slightly smaller) of the exponent observed in the scaling relation of figures ^ and 
H which, respectively, predict values of 0.85 ± 0.09 and 0.86 ± 0.10 for x. 
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Typically, the observed scaling dependencies in these types of problems result 
from the local stress balance. A scaling analysis where the viscous stresses 
of the upper and lower fluids balance the stress arising from the interfacial 
curvature predicts linear scaling dependencies and conical profile shapes. The 
non-linearity of the observed dependencies and the lack of dependence of the 
similarity solution on i/j indicate that a different stress balance may govern the 
flows (e.g. only viscous stress due to upper fluid balances stress due to the 
interface curvature). Furthermore, non-local effects coupling themselves into 
the solution could account for the slight differences between the values of the 
exponents /3 and x. Further evidence for the existance of these coupling effects 
is presented in section ^. 




Figure 9: The scaled hump profiles for system 6. The lower inset shows eight 
profiles taken from the S = 0.830 data set. The main figure shows the same 
profiles after scaling. The solid line corresponds to a power law of the form 
In the upper inset we compare the universal curves for the S — 0.830 cm, 
0.613 cm, 0.508 cm, 0.381 cm, 0.255 cm data sets. This figure is reproduced 



from Cohen & Nagel (?) 



In both figures, the upper right inset shows a comparison of the similarity 
solutions for five different tube heights. The profiles corresponding to the dif- 
ferent tube heights all display the same power-law dependence. Within error, 
^ (taken from figures ^ and |8|) can be used to scale the radial components of 
these profiles and obtain good collapse. Recall that the normalization prefac- 
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Figure 10: The scaled hump profiles for system 5. The lower inset shows eight 
profiles taken from the S = 0.667 data set. The main figure shows the same 
profiles after scaling. The solid line corresponds to a power law of the form 
R^-^^. In the upper inset we compare the universal curves for the S = 0.984 
cm, 0.921 cm, 0.889 cm, 0.667 cm, 0.445 cm data sets. 
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tors n{S) decrease roughly as exp"^-^"^. Here, this decrease is correlated with 
the observation that the profiles become shallower at larger S. The points of 
deviations for the S = 0.255 cm and 0.381 cm profiles in the large vi system, 
and the S = 0.445 cm profile in the small vi system, mark the transition from 
the similarity regime to the matching regime beyond which the profiles become 
asymptotically flat. At large enough radii all of the scaled profiles display these 
deviations. 

As a final check, figure [l| shows a comparison of the similarity profile for 
S = 0.831 cm in the large fi system and the similarity profile for S = 0.667 
cm in the small vi system. An error analysis calculatiorj^ for the data points 
located between — 5 < i? < 5 (where the residuals are centered around zero - 
see figure |l]) shows that ~ This value indicates that the differences 
between the two similarity curves are on the same order as the experimental 
uncertainty which is taken to be a quarter of a scaled pixel. The excellent col- 
lapse verifies that I'l and the viscosity ratio do not affect the similarity solution 
and emphasizes that the scaling behavior observed in these systems is robust. 



5 Mapping Out The Transition Location 

Having shown that the viscosity ratio does not affect the scaling relations and 
self-similarity characterizing the detailed structure of the transition, I return to 
the question of whether Su is affected by such a change in parameters. There is a 
large body of work dating back to the late 1940's (see, for example. Rouse (19561 ); 



Gariel (1949| ); praya (1949| ); [Muskat (19"72| )) which addresses the mapping out 



of the selective withdrawal transition location. A large portion of this work 
focuses on the problem of extracting crude oil deposits without withdrawing 
any of the water which is often trapped beneath the oil. Nearly all of the 
experimental studies found in the literature which address these large scale 
extraction problems assume that the flows have a high Reynolds number or, 
equivalently, that the viscous stresses are negligible. Therefore, low viscosity 
fluids were used to model the flows. One exception arises in the modeling by 



Blake &: Ivey (1986 ) of magma layer mixing during volcanic eruptions. However, 



as described earlier, miscible fluids were used in those investigations. Since, as 



Lister (1989 ) showed, the absence of surface tension allows for the withdrawal 
of both fluids at any withdrawal rate, these particular studies never addressed 
the actual selective withdrawal transition. With the advent of new technologies 

^ Since the data points for both curves were not aligned, wa^s calculated in the following 
way: First a local linear fit was used to interpolate the value of curvei between the data 
points. Second, the minimum distance between the interpolated curve and the value of the 
points in curve2, Si^2 was calculated. The value for the entire curve is defined as: = 

'i I 'i \/ where N is the number of points compared, crix and aiy are 

the experimental errors associated with the x and y points in curvei, cr2x and a2y are the 
experimental errors associated with the x and y points in curve2 , n is equal to the number of 
points used in making the linear fit to curvei, and the sum is taken over the index i which 
labels the points in curve2. 
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Figure 11: Comparison of the similarity profiles for systems 5 and 6. The 
closed symbols correspond to the S = 0.830 cm data for system 6 while the 
open symbols correspond to the S = 0.667 cm data for system 5. Figure pl] b 
shows a plot of residual quantity, Si for all of the data points. Here, Si is the 
minimum distance between each point on the system 6 curve and the value 
of a local linear fit to the system 5 curve at the same R value. The region 
between the vertical dotted lines corresponds to the similarity regime where 
the residuals become centered around zero. The calculation was performed 
within this similarity regime. 
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which use the selective withdrawal geometry in conjunction with low Reynolds 
number flows (see, for example, Cohen et al. (2001); Ganan-Calvo (199^ )), it 
has become increasingly important to determine the effect of the viscous stresses 
on the transition. 

Figure |l^ shows a plot of Su versus Q for seven pairs of fluids. Table |l| lists 
the fluids used in making these measurements along with the measured values 
of the interfacial tension, 7, the fluid densities, pu and pi, the density mismatch, 
A/9, the fluid viscosities, and vi^ and the viscosity ratio vi/vu- 




10"^ 10"^ 10° 10^ 
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Figure 12: Plots of the transition tube height Su as a function of Q for the seven 
systems listed in table |l|. Table |^ lists the symbol used in plotting the Su curve 
for each system. For each system, it is observed that Su oc where the power 
a (listed in table |l|) ranges between 0.30 and 0.45. The Su curves cluster into 
three groups which correspond to the three different upper fluids used in the 
experiments. Each cluster is labeled with the corresponding upper fluid. 

There are a few obvious trends. First, all of the data sets show that Su oc Q" 
where the power a ranges between 0.30 and 0.45. The last row in Table | 
shows the value of the exponent for the power-law used to fit the Su curves 
for each system. Second, for a given value of Q, an increase in the upper 
fluid viscosity, increases Su- The three distinct clusters of curves in figure ^ 
correspond to three distinct upper fluid viscosities. As figure |l2| shows, a five 
fold increase in Uu can increase Su by a factor of two. Third, for a given value of 
Q, even a thousand fold increase in the lower fluid viscosity docs not significantly 
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affect Su- These last two observations indicate that it is the upper fluid viscosity 
rather than the viscosity ratio which affects the transition location within the 
S vs. Q parameter space. 

The remaining trends in the data are weak. In order to ascertain the effects of 
the surface tension, 7, and density mismatch, Ap on the transition straw height, 
these parameters must be varied by large amounts. Unfortunately, reducing 7 
below about 20 dynes/cm causes mixing of the fluid layers and results in a 
diffuse interface under high shear rates. Therefore it is very difficult to vary 7 
in the experiments by more than a factor of two. Another difficulty arises from 
an inability to decouple the changes in these parameters for different systems. 
For example, when comparing systems 1 and 4, an order of magnitude increase 
in Ap is coupled with a factor of two decrease in 7. Since both Ap and 7 
act to stabilize the interface, the effects due to these variations may cancel. A 
further complication arises due to the effect of surfactant concentration on the 
interface. As will be discussed in Appendix A, the accumulation of surfactants 
on the interface tend to shift Su by an amount which is comparable to the slight 
shifts seen in top four curves (which form the uppermost cluster) in figure 
More careful measurements of Su are necessary in order to resolve these weaker 
trends. 

It is instructive to compare these results with currently available theoretical 
predictions for the transition which take into account both the viscous stresses 
and the surface tension. While there are no current numerical studies which 
have been tailored to address this exact problem, the closest approximation to 
such a treatment can be found in the Lister (1989| ) paper on selective with- 
drawal for zero Reynolds number flow. While his simulations were performed 
for equal viscosity fluids, the fact that the lower fluid viscosity is irrelevant to 
the transition location can be used to compare the simulations with the exper- 
iments for which the fluid viscosities are unequal. Since those simulations were 
designed to address large scale magma flows rather than the small scale fluid 
flows found in the experiments, care must be taken in making a comparison 
of the results. For example, the simulations were designed to model a system 
with a point sink which is located many capillary lengths (defined as \J -^^) 

away from the interface. However in the experiments, S is comparable to the 
capillary length. Thus, the experiments are performed at much smaller values of 
Q and Su than the simulations making it impossible to quantitatively compare 
the results. Nevertheless, since it is expected that Su is a smooth function of 
Q, one can check that the two data sets are consistent with this expectation. 

Another issue which needs to be addressed is the determination of the 
Reynolds number for the experimental flows. There are a variety of ways in 
which the Reynolds number can be deflned. The greatest change in the velocity 
of the upper fluid occurs along the line connecting the tip of the hump and the 
withdrawal tube. Since the fluid at the tip of the hump is stationary the magni- 
tude of this change in velocity is ^(^2)^ ■ '^'^'^ distance of the tube oriflce to the 
interface is the largest length scale affecting the flowsQ. Therefore, the quantity 

*For this range of flow rates and density mismatch values, the capiUary length scale is a 
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R^max = ■nipfiy^v used as an upper bound for the Reynolds number. 

In figure |l2[ the upper most cluster of Su curves has 0.005 < Rcmax < 6, the 
middle cluster of curves has 0.03 < Rcmax < 50, and the lowest curve has 
0.05 < Re„iax < 350. These estimates imply that a comparison with results 
from a zero Reynolds number simulation must be treated with caution. For 
example, the effects of the fluid's inertia could manifest themselves as a change 
in a, the value of the exponent for the power-law fits in figure |l2|^. Such an ex- 
planation could account for the decrease in a which occurs when i^u is decreased 
from 10 St to 1.9 St. However, since an equally large change in Reynolds num- 
ber can occur along the Su versus Q curves it is curious that kinks which may 
be indicative of such a change are not observed. 

In figure |l^ the Su vs. Q curves are compared with the predictions of Lister 
for systems with i^u — 10 St, i^u = 1-9 St, and Vu = 0.57 St corresponding 
to systems 4, 6 and 7 in table |l|. The solid and open symbols denote data 
from the experiments and simulations respectively. The dashed lines represent 
power-law fits to the experimental results. These fits are projected into the 
regime for which the simulations have a prediction for the value of Sw Since 
Lister's simulations were performed using dimensionless variables, the proper 
values for the fluid parameters were used to determine the values of S and Q. 
Furthermore, the simulation flow rate is reduced by a factor of two in order to 
make a first order correction for the fact that the simulations use a point sink 
rather than a tube to withdraw the fluids. 

There are two features of the comparison which need be addressed. First, 
the simulations curves for the f„ = 10 St system fall below the projected power- 
laws. In Section H evidence was provided for the lack of dependence of Su on the 
straw diameter D for the range of flow rates used in the experiments. Therefore 
it is unclear whether geometric factors alone could account for the large shift 
which is necessary for the collapse of the numerical results onto the projected 
experimental power-laws. Second, it is disturbing that the power-law exponents 
for the systems which use less viscous (rather than more viscous) upper fluids 
agree more favorably with the zero Reynolds number simulation results. As 
stated above, the shift in the power-law exponent for the experimental results 
indicates that there may be more than one way to balance the stresses acting 
on the interface. The simulation results only address the regime where surface 
tension and buoyancy effects play an equal role in stabilizing the interface and 
therefore do not observe the change in the S vs. Q power-law exponent^ 
Nevertheless, the fact that the simulations predict a power law exponent of one 
third and the fact that the experimentally measured exponents are centered 
around this value deserves notice. 

Finally, it is often possible to use the scaling relations observed in systems 
approaching a singularity to help identify which stresses are balancing each 
other. Such an identification for the different fluid systems could in turn help 

little smaller than S. However, it is conceivable that for a set of density matched fluids or for 
measurements at very small S the capillary length could become larger than S. 

^W. W. Zhang private communication. 

''W. W. Zhang private communication. 
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Figure 13: Comparison of the experimentally measured Su curves for systems 4, 
6, and 7 in Table |^ with those predicted by the simulations of Lister (1989). The 
closed symbols depict the experimental measurements while the corresponding 
open symbols depict the simulation predictions for a system with the same fluid 
parameters. The dashed lines are power-law fits to the experimental data. 
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explain the observed changes in the exponent a. 

6 Conclusions 

In this paper, it has been shown that the transition straw height for a given 
flow rate, Su, changes as: Su oc Q", where a ranges between about 0.45 and 
0.30. This is the first instance in which experiments comparing Su for different 
systems of immiscible fluids with flows in the low Reynolds number regime have 
been performed. The data indicates that for the selective withdrawal system, 
the upper fluid viscosity rather than the viscosity ratio determines the value of 
Su- It has been shown, that when the withdrawal tube is sufficiently close to 
the interface, the transition is hysteretic: the straw height at which the spout 
collapses back into a hump is different from S'„. At larger straw heights, the 
difference between Su and the spout to hump transition straw height, AS, decays 
exponentially. This effect has been linked to the tube diameter which at low S, 
sets a length scale for the final mean radius of curvature l/nu- However, the 
transition remains discontinuous and Ku remains finite even at large S. 

I have performed a detailed comparison of the scaling relations for two sys- 
tems which have the same upper fluid viscosity but have a lower fluid viscosity 
which is different by a factor of two hundred. The surface tension for the two 
systems was different by about 13% while the density mismatch was different 
by a factor of three (See systems 5 and 6 in table |^). It is observed that up until 
the cutoff, the hump proflles behave as though they are approaching a singular 
solution where, at the flow rate Qc the hump height would be equal to and 
the mean curvature, k, would diverge. The quantity ^"^^""^ has been shown to 
scale as {^)'^ where /3 = -0.85±0.09 for the high i/, system and (3 = -0.86±0.10 
for the low vi system. These scaling relation were used to collapse the hump 
profiles for different flow rates and straw heights near the transition onto a sin- 
gle universal curve. The region of the similarity profiles located beyond the 
parabolic tip can be fit with a power law which has an exponent x = 0.72 for 
both systems. The results show that both the scaling exponents and the shape 
of the similarity solution are independent of i^i and the viscosity ratio. In fact, 
a direct comparison of the similarity solution for both fluid systems indicates 
that, within error, the curves are identical. 

As described in the introduction, in the 2-D analogue to the selective with- 
drawal problem the lower fluid viscosity plays a fundamental role in determining 
how close the system can approach the cusp singularity. In the 3-D selective 
withdrawal system, the asymptotic value of Ku shows little or no dependence on 
ui. Therefore, the effects of the lower fluid viscosity must enter as higher order 
terms which are undetected by the experiments. Since k„ is independent of the 
flow rate at large S, an increase in the Capillary number, Vu^/j can also be 
ruled out as a method of getting the system closer to the singularity. 

The saturation value of the mean curvature, Kusat does show some depen- 
dence on the upper fluid viscosity Uu- For example, systems 5 and 7 in table || 
show that when 7 and Ap are kept constant, decreasing ^'„ from 2.0 St to 0.57 
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St, decreases Kusat by a factor of six. However, as systems 2 and 5 in the table 
indicate, this trend is not uniform. Shallow hump profiles are also observed for 
nearly inviscid systems which use air as the upper fluid and water as the lower 
fluid. Therefore, it is likely that one of the ways in which nw„ can affect Kusat 
is by shifting the balance between the viscous and inertial terms in the gov- 
erning Navier-Stokes equations. For the remaining systems which have a lower 
Reynolds number characterizing the flows, even a factor of five increase in nu^ 
does not significantly affect Kusat- 

Since at low enough Reynolds numbers the fluid viscosities do not affect 
Kusat 1 tbe only remaining parameters which may act to set the length scale 
for the cutoff are the density mismatch Ap and surface tension 7. With the 
exception of systems 1 and 6, the large error in the data makes it difficult to 
distinguish between the Kusat values for the different systems. Also, the trends 
in the data are too weak to allow for a determination of the effects of 7 on Kusat ■ 
However, a comparison of system 1 with systems 2, 3, and 4 and a comparison 
of system 6 with system 5 indicates that an increase in Ap results in a higher 
value for Kusat- Note that for both comparisons Vu is kept constant. This 
result is intriguing since the similarity treatment for the profiles is localized to 
the hump tip whereas the effects of the density mismatch should only affect 



the system on a length scale comparable to the capillary length y These 

observations suggest that the matching region which connects the profile near 
the tip of the hump to the flat interface at large radii is responsible for setting 
the length scale for the cutoff f\. Such an effect could also explain the slight 
difference between the observed value of the power-law exponent describing the 
similarity solution(x) and the value predicted by the scaling relations(/?). The 
importance of identifying which parameters determine the length scale for the 
cutoff warrants a more careful investigation of the kusat dependence on Ap, 7 
and the local interfacial boundary conditions. 

The robustness of the similarity analysis shows that singularities can be used 
to organize the study and classification of the steady state hump profiles near 
the selective withdrawal transition. In particular, the discontinuous nature of 
the transition, marked by the cut-off curvature Ku, coupled with the display 
of scaling behavior suggests a transition structure which is remarkably similar 
to that of wcakly-first-order thermodynamic transitions. Whether this analogy 
hints at some deeper relationship between the classification schemes for weakly 
first order thermodynamic transitions and those for these types of topological 
transitions remains to be shown. 
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editing this manuscript. This research is supported by the University of Chicago 

'^J. Eggers, W. W. Zhang, S. R. Nagel, J. Wyman, and H. A. Stone (private communica- 
tion) . 
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A Surfactant effects 

Surfactants are chimeric molecules consisting of hydrophobic and hydrophilic 
sections which tend to aggregate at fluid interfaces^. Even under very controlled 
conditions, it is difficult to keep fluid systems free of surfactants for long periods 
of time. Since, surfactants can significantly affect the interfacial tension and 
surface flows, part of understanding the details of any fluid interface problem 
entails isolating and accounting for the effects of surfactants. 

In the present studies there are two types of effects which are thought to 
result from the presence of surfactants. First, for a given flow rate, it is observed 
that over a period of several days, 5'„ increases by a small amount indicating that 
the surface tension in the system may have been uniformly reduced. Second, 
at the transition, the system oscillates between the hump state and the spout 
state over a period of about minuten in a manner similar to that observed by 



de Bruijn (1993) in tip streaming of drops under shear flow. In the tip streaming 
experiments, a straining flow stretches a drop of fluid so that the two tips of the 
drop (each of which is analogous to the hump tip) have a high curvature. At 
sufficiently high shear rates the drop can enter the streaming state (analogous 
to the spout state) where a small jet of fluid eminates from each of the drop 
tips. It is observed that near the transition the drops oscillate between the 
stretched drop state and the streaming state. This oscillatory behavior has been 
linked to the accumulation of surfactants near the stagnation point located at 
the drop tip. This accumulation lowers the local surface tension and causes a 
transition to the streaming state. The jet "sweeps" the interface, reduces the 
local concentration of surfactants, and subsequently collapses back into the drop 
state at which point the cycle begins again. While the local boundary conditions 
for this problem are different from those of the selective withdrawal problem, 
the flow patterns are remarkably similar indicating that the oscillatory behavior 
observed in both systems may be correlated. Below, I explain how both of these 
effects influence the results presented in the main body of this paper. 

Figure ^ shows plots of Su vs. Q for three data sets corresponding to 
measurements performed on system 6. The filled down-triangles correspond 
to measurements taken just after the fiuid interface was cleanedFn The open 



down-triangles, correspond to measurements of Su taken after the surfactant 

*For some oil-oil (for example Polybutadiene and PDMS) interfaces the random thermal 
energy is able to keep the surfactants off the interface. However these systems usually have a 
very low surface tension and therefore develop diffuse interfaces under high shear rates. 

^Note that the time scale for this effect is much larger than the time scale for the pump 
induced noise in the withdrawal rate which can produce a similar effect when the system is 
very close to the transition. 

^"When the system is in the spout state, surfactants, which have accumulated on the inter- 
face are "swept" into the straw and deposited in the waste container. By leaving the system in 
the spout state for long periods of time, the surfactant concentration over the entire interface 
can be significantly reduced. 
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concentration at the interface was allowed to equilibrate for a period of a week. 
A uniform 20% increase in Su, or, equivalently, a 50% decrease in the transition 
flow rate is observed. Upon cleaning the interface once again, the Su data points 
return to their original value (open diamonds). Note that the observed shift in 
Su resulting from the differing surfactant concentration is quite small compared 
with the shift observed when i/u is increased by a factor of five. Nevertheless, 
the fact that the transition flow rate can change by nearly 50% depending on 
the surfactant concentration is noteworthy. Finally, in figure a comparison 
between the ^'=^^"°°' vs. ^ curves for the equilibrated system (squares) and 
the clean system (open circles) indicates that the shift in Su does not affect the 
scaling relations. 
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Figure 14: Plots of the transition tube height Su as a function of Q for sys- 
tem 6 in Table |l} The filled down-triangles correspond to measurements taken 
just after the interface was cleaned. The open down-triangles correspond to 
measurements taken after the surfactant concentration at the interface was al- 
lowed to equilibrate over a period of a week. The open diamonds correspond to 
measurements taken after the system was cleaned once again. 

While the increase in the Su values is most likely due to a uniform reduction 
in 7 over the entire interface, the dynamic nature of the observed hump to 
spout oscillations implies that (as with tip streaming) the flows may cause local 
variations in the surfactant concentration. In this picture, when the system is 
in the hump state, surfactants are dragged towards the hump tip by the surface 
flows. If the system is sufficiently close to the transition, this accumulation, 
which lowers the surface tension locally, causes the hump to increase its height 
and curvature and ultimately drives the system into the spout state. Once in 
the spout state, the local surfactant concentration reduces over time and the 
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Figure 15: Plots of 



''-'^^ vs. ^ for system 6 in Table |l]. The squares corre- 
spond to measrements taken after the surfactant concentration at the interface 
was allowed to equilibrate over a period of a week. The open circles correspond 
to measurements taken after the system had been cleaned. 



system eventually decays back into the hump state. The detailed characteristics 
of this cyclical behavior depend on the surfactant concentration over the entire 
interface. 

When the entire interface has been cleaned, it takes a long time for the sur- 
factants to accumulate in concentrations large enough to affect the steady state. 
Therefore, the increase in the hump curvature and height is slow. Furthermore, 
once the interface forms a spout, it takes a very short amount of time for the 
local concentration of surfactants to be reduced. Consequently, systems that 
have been cleaned spend a very small fraction of the hump to spout oscillation 
cycle in the spout state. Since many of the measurements discussed in the paper 
were performed in the vicinity of the transition, care needed to be taken when 
measuring the hump height and curvature. In order to understand how hmax 
and K change with time, the following experiment was performed. First, the 
system was placed in the hump state. Then, the flow rate was increased mo- 
mentarily so that the interface formed a spout for a short period of time. Just 
after spout collapse, the time dependence of the hump curvature and height was 
measured. Figures and show the results of measurements performed 
for the same system, with the same value of S, but at different values of Q. 
The initial decay, which occurs over a time scale of about five seconds, gives 
some measure of the relaxation time for the flows in these systems. Following 
this initial decay, the height and curvature values plateau. These plateaus are 
associated with a regime where the local surfactant concentration is too low to 
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affect the shape of the interface. When the system is near the transition so that 
^"q'^ « 0.004, the plateau regime is short hved (about fifteen seconds) and is 
followed by a regime in which both the height and curvature increase their val- 
ues (figure p^). However, when the experiments are performed at larger values 
of '^"g'^ the plateau regime is longer. Figure shows that when ^"q"^ ~ 0.01 
the plateau regime lasts for over twenty minutes. 

When collecting data for the scaling relations, it is important to determine 
the value of the hump height and mean curvature in the plateau regime. Also, 
a larger amount of error must be assigned to data points taken from plateau 
regimes that are short lived. Note however that since the plateau regime length- 
ens quite rapidly as ^ is increased, for the cleaned interface, these precau- 
tions only apply to the one or two data points in the scaling relations, which 
are closest to the transition. 

While all of the scaling data presented in the main body of this paper was 
taken for clean interfaces, it is useful to understand the changes in the transition 
structure which occur when the surfactant concentration has been allowed to 
equilibrate over a period of days. For equilibrated interfaces the increase in 
the hump height and curvature occurs over a shorter amount of time so that 
the plateau regimes shorten. Therefore, greater caution needs to be taken in 
the measurements of the scaling relations. However as figure ^ shows, when 
such precautions are taken no noticeable changes in the scaling relations are 
observed. 

Another characteristic of the equilibrated systems, is that (unlike the cleaned 
systems) they can spend a large fraction of the hump to spout oscillation cycle in 
the spout state. During the time which the system spends in the spout state, it 
is possible (even at high Q or S) to quickly raise the straw and cause the spout 
to collapse back into the hump state|^. The amount AS' necessary to cause 
the collapse can be measured. I find that the value of AS varies depending on 
conditions. For example, if the AS measurement is performed just after the 
local concentration of surfactants has been reduced, only a small amount of 
hysteresis is observed. On the other hand if the experiment is performed after 
the concentration of surfactants has been allowed to build up, AS can become 
as large as 80 ^m. Note that even this large amount of hysteresis is orders of 
magnitude smaller than the typical value of the straw height and therefore does 
not significantly affect the Su vs. Q curves for the equilibrated system. 

In conclusion, it is observed that allowing the surfactant concentration to 
equilibrate over a period of days has a relatively small effect (20% increase) on 
the Su measurements (figure ^4|) and, when the propper precautions are taken, 
no effect on the scaling relations (figure p^ . 

Since the spout drains the lower fluid, the straw height S increases with time. Therefore, 
it is necessary to check that the amount of hysteresis is actuaUy larger than the increase in S 
which occurs during the time it takes to perform the measurement. 
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Figure 16: Plots of the time dependence for hmax (squares) and k (open circles) 
just after spout collapse for system 3 in Table |l|. The figures correspond to 
measurements taken at constant S but at different Qs. In both figures, to 
corresponds to the time at which the spout collapses into a hump. Figure 16 1 



shows measurements for a system which is very close to the transition with 
'^"q^ ~ 0.004. Figure shows measurements for a system which slightly 

farther away from the transition with ^"q^ ~ 0.01. 
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